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ABSTRACT 

A few Galactic classical Cepheids were observed in the programmes of space missions 
as Coriolis, MOST and Kepler. An appealing opportunity was to detect additional non- 
radial modes, thus opening the possibility to perform asteroseismic studies and making 
the pulsational content of Galactic Gepheids more similar to that of Magellanic Clouds 
ones. However, only hints of cycle-to-cycle variations were found, without any strict 
periodicity. In this context the potential of the CoRoT exoplanetary data base was 
not fully exploited despite the wide area covered on the Galactic plane. Therefore, we 
investigated all the candidate Cepheids pointed out by the automatic classification of 
the CoRoT curves. At the end we could identify seven bona-fide Cepheids. The light 
curves were investigated to remove some instrumental effects. The frequency analy¬ 
sis was particularly delicate since these small effects can be enhanced by the large 
amplitude, resulting in the presence of significant, but spurious, peaks in the power 
spectrum. Indeed, the careful evaluation of a very attracting peak in the spectra of 
CoRoT 0102618121 allowed us to certify its spurious origin. Once that the instru¬ 
mental effects were properly removed, no additional mode was detected. On the other 
hand, cycle-to-cycle variations of the Fourier parameters were observed, but very small 
and always within ±3cr. Among the seven Cepheids, there are two Pop. I first-overtone 
pulsators, four Pop. I fundamental mode pulsators, and one Pop. II star. The CoRoT 
colours allowed us to measure that times of maximum brightness occur a little earlier 
(about 0.01 period) at short wavelengths than at long ones. 

Key words: methods: data analysis - stars: individual: V376 Mon - stars: individual: 
DP Mon - stars: interiors - stars: oscillations - stars: variables: Cepheids 


1 INTRODUCTION 

The use of Cepheids as stellar candles requires a deep knowl¬ 
edge of their pulsational characteristics. The separation be¬ 
tween fundamental (F-mode) and first-overtone (lO-mode) 
radial pulsators (Antonello, Poretti & Reduzzi 1990) is still 
the most relevant distinction to be taken into account to cal¬ 
ibrate the period-luminosity relations. About the regularity 
of their variations, ground-based photometry of V473 Lyr 
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(Burki, Mayor & Benz 1982) and Polaris (Evans, Sasselov 
& Short 2002) showed continuous changes in amplitude and 
phase. Molnar & Szabados (2014) established that V473 Lyr 
pulsates in the second radial overtone and shows strong am¬ 
plitude and phase modulations on two long periods (1205±3 
and 5300±150 d). The slow changes in the light curve of Po¬ 
laris on a time baseline of 4.5 yr were precisely recorded by 
space photometry using the Coriolis spacecraft. Spreckley 
& Stevens (2008) connected them to a quite fast evolution¬ 
ary phase. However, Bruntt et al. (2008) proposed a differ¬ 
ent interpretation pointing out how the amplitude change is 
cyclic rather than monotonic and also suggesting an excess 
of power in the 2-6 d range due to the granulation. 

The situation was expected to change with the advent 
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of high-precision space mission. A detailed analysis of the 
MOST observations of SZ Tan and RT Aur revealed vari¬ 
ations from cycle to cycle as a function of the pulsation 
phase (Evans et al. 2015). The stars were both monitored 
for less than seven consecutive cycles: the lO-Cepheid SZ 
Tan showed a pulsation less stable than the E-Cepheid RT 
Aur. The historical analysis of the times of majcimum bright¬ 
ness proved that both stars experienced period variations on 
a time-scale of decades and that those of the overtone mode 
resulted to be more erratic than the F — mode ones. 

In the case of Kepler, there was only one Cepheid in the 
field-of-view, i.e., V1154 Cyg. Derekas et al. (2012) provided 
a very detailed analysis of the data acquired in the first 600 d 
(i.e., about 120 cycles) of the mission, detecting significant 
cycle-to-cycle fluctuations. A very slight correlation was also 
found between the Fourier parameters and the 0-C values, 
i.e., the difference between the observed and calculated times 
of maximum brightness (Tmax). This suggested that the O- 
C variations (up to 30 min) might be due to instabilities in 
the light curve shape. On the basis of these observational 
results, Neilson & Ignace (2014) suggested that convection 
and hot spots can explain the observed period jitter, thus 
supporting the Bruntt et al. (2008) interpretation of the Po¬ 
laris changes. The clues about the changes in the light curve 
of VI154 Cyg were strengtheued by the analysis of new Ke¬ 
pler measurements (up to Q17 quarter): the analysis of the 
Fourier parameters supplied a robust determination of a pe¬ 
riodicity of 158.6 d (Kanev, Savanov & Sackhov 2015). 

About ground-based observations, Anderson (2014) dis¬ 
covered significant modulations in the radial velocity curves 
of some Cepheids. These modulations appear to be different 
from those reported on V1154 Cyg, since they manifest as 
shape and amplitude variations on time-scales of years for 
short-period Cepheids (like VI154 Cyg) and cycle-to-cycle 
variations in the long-period Cepheids. 

A common result among all these investigations is that 
no nonradial mode was identified in the light or radial ve¬ 
locity curves of Galactic Cepheids. Hence the counterparts 
of the subclasses detected in the Large Magellanic Cloud 
(LMC; for a review see Moskalik 2014) are still missing. 
However, an important player has not yet contributed, i.e., 
the CoRoT data base of the light curves collected in the 
exoplanetary fields. Such data base already proved to be a 
real treasure allowing the discovery of a unique triple-mode 
Cepheid in the outskirts of the Milky Way (Poretti, Baglin, 
& Weiss 2014). In this paper we describe the new elements 
brought by the CoRoT photometry to the hot topic of the 
stability of the Cepheid light curves. 


2 COROT DATA AND SPECTROSCOPIC 
FOLLOW-UP 

The first new Cepheid discovered by CoRoT was 2MASS 
J06415168-012005feCoRoT 0102618121, a 12“'-mag star lo¬ 
cated at a = 06*'dP” 5F685, 5 = -01° 20'05'.'960 (J2000). 
It was observed during the first Long Run in the Anticentre 
direction (LRaOl) and then in the LRa06 four years later 
(Table 1). The classification of the type of variability was 
a little troublesome. The star was assigned to the RR Lyr 
Working Group, but the preliminary analysis immediately 
pointed out that the star was not a RR Lyr variable since 


the period was too long. The Cepheid variability was much 
more probable, but some doubts still remained, based on 
the classification as MIV type reported in the header of the 
FITS file of the LRaOl data. Indeed, the “CoRoT Variability 
Classifier” automated supervised method (CVC; Debosscher 
et al. 2009) suggested an active, rotational variable as first 
hypothesis and a classical Cepheid as a second one. The 
spectral classification of CoRoT 0102618121 was revised as 
A5V in the current version of the Exo-Dat (Deleuil et al. 
2009) catalogue and consequently reported in the header of 
the more recent LRa06 data. 

While the analysis of CoRoT 0102618121 was in 
progress, we proposed known Cepheids as CoRoT tar¬ 
gets by answering the Announcements of Opportu¬ 
nity regularly issued before any run. In such a way 
V376 Mon=CoRoT 0300003751 and DP Mon=CoRoT 
0221644967 were monitored in the LRa02 and SRa02 runs, 
respectively (Table 1). However, it should be expected that 
other Cepheids would be discovered thanks to the deep sur¬ 
vey performed by CoRoT on the Galactic plane. 

To search for these new Cepheids, we first scrutinized 
all the new variables discovered with the Berlin Exoplanet 
Search Telescope (BEST) that have been classified as can¬ 
didate Cepheids (Karoff et al. 2007; Kabath et al. 2007, 
2008, 2009a,b; Fruth et al. 2012) and actually observed with 
CoRoT. We noticed that the variability of V376 Mon was 
recorded (Kabath et al. 2009b) with the correct period, i.e., 
1.652 d (Berdnikov et al. 2014), not the 0.623 d reported in 
the discovery paper (Wachmann 1964). Unfortunately none 
of other BEST candidates were confirmed as a new Cepheid, 
since they turned out to be rotational variables. The misclas- 
sification is probably due to the too relaxed criterion used to 
classify a new variable as Cepheid, i.e., a periodic variability 
with an amplitude larger than 0.10 mag. Indeed, such an 
amplitude threshold can be reached by active spotted stars. 

As a second step in the data mining of the CoRoT 
data base, we selected all the variables for which the CVC 
method proposed a classification as Cepheid, independently 
from rank and probability. The full light curves of about 
250 candidates thus selected were inspected by eye. Eclips¬ 
ing variables, RR Lyr and high-amplitude 5 Set stars were 
found and not considered furthermore. We remained with a 
sample where the dichotomy was again between Cepheid 
and rotational variables. The stars characterized by very 
large cycle-to-cycle variations were immediately considered 
as rotational variables. The stars showing more stable light 
curves were aualysed in frequency and least-squares solu¬ 
tions of the mean light curves were obtained. We note that 
in this process the small amplitude was not considered as a 
rule of thumb since it could be reduced by an unresolved, 
blended star. On the other hand, a large amplitude (say, 
greater than 0.20 mag) was recognized as the signature of a 
pulsational mechanism. Most of the stars were considered as 
rotational variables after that an ensemble of several param¬ 
eters (period, amplitude, stability of the light-curve shape, 
Fourier parameters, ...) was taken into account. However, we 
cannot rule out that a few of them are actually Cepheids, 
whose regular curve is perturbed by that of a blended irreg¬ 
ular variable. We could identify some borderline cases (e.g. 
CoRoT 0102635101, 0102640143, 0605539616, 0223971984, 
0102578739, 0110827792, and, overall, 0223971984) where 
the distinction between activity and pulsation was very difB- 
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Figure 1. Frequency analysis of CoRoT N2 data of 
CoRoT 0102618121. Upper part, LRaOl data: example of light 
curve (top panel), residuals after prewhitening of the main peri¬ 
odicity from the global dataset (middle panel), behaviour of the 
mean flux all along LRaOl (bottom panel). Lower part-, same as 
in the upper part, but for LRa06 data. 

cult. New spectroscopic observations are necessary to estab¬ 
lish the nature of these problematic stars in an unambiguous 
way. 

During the screening procedure we discovered the 
unique Galactic triple-mode Cepheid CoRoT 0223989566 
(Poretti et al. 2014) and at the end we could count seven 
bona-fide Cepheids. In addition to CoRoT 0102618121, 
V376 Mon, and DP Mon we identified four new 
Cepheids: CoRoT 0221640090, 0315221415, 0211626074, and 
0659466327 (Table 1). The CoRoT photometry of V376 Mon 
and of other Cepheids was analysed when the interpre¬ 
tation of the CoRoT 0102618121 data was almost fin¬ 
ished. Therefore, we describe the procedure as applied to 
CoRoT 0102618121 first and then how it has been extended 
to other Cepheids later. 

3 THE ANALYSIS OF CoRoT 0102618121 

DATA 

3.1 Spectroscopic follow-up 

To solve the ambiguous spectral classification of 
CoRoT 0102618121 and to validate our classification 
criteria, we obtained three high-resolution spectra on 
2012 January 10 and 11 with the HARPS spectrograph 
installed at the 3.6-m telescope of the ESO-La Silla Ob¬ 
servatory. The spectra were taken in the high-efficiency 
EGGS configuration (R=80,000) and they have signal-to- 
noise ratios (S/N) between 50 and 55. We reduced and 
normalized the spectra with a semi-automated pipeline 
developed at the Brera Observatory (Rainer 2003). We 


then computed the mean line profiles for each spectrum 
using the LSD software (Donati et al. 1997) and we used 
them to estimate the radial (from 47.9 to 53.6 kms“^) and 
projected rotational (from 9.64 to 11.85 kms“^) velocities 
by fitting the profiles with gaussians. The reduced and 
normalized spectra were shifted by their radial velocities 
and averaged in order to lessen the effect of the pulsation. 
The spectral synthesis software SME version 4.2.3 (Valenti 
& Piskunov 1996) was used to estimate the physical 
parameters. The spectral regions 5160-5190 A (including 
the gravity sensitive Mg-triplet) and 4830-4900 A (the H/3 
region, sensitive to the temperature) were selected. We left 
temperature (Teff), gravity (logg, cgs units), metallicity 
([Fe/H]), micro- and macroturbulence as free parameters, 
while usini = 11 kms“^ was imposed. Following Valenti & 
Fischer (2005), we also left the abundances of Na, Si, Ti, 
Fe, and Ni as free parameters. We performed an iterative 
process on the MARCS stellar models (Gustafsson et al. 
2008): we started running the SME software on the Mg 
region, then we switched to the H/3 region keeping the 
logg previously found as a fixed parameter. After that 
we went back to the Mg region keeping fixed the 
value found in the region, and so on until the results 
converged. The resulting parameters are Tf,g= 6087T44K, 
logg=1.53±0.15 dex, [Fe/H]=-0.34 ± 0.04 dex. They 
confirm that CoRoT 0102618121 is a Pop. I supergiant and 
not a main-sequence star: its classification as a Cepheid 
variable was definitely established. 

3.2 CoRoT data 

The original timeseries of CoRoT 0102618121 consist of 
348 702 and 206 528 data points, respectively. In both cases 
the exposure time was set to 32 s (corresponding to the 
CoRoT short cadence) and the flux was recorded in three dif¬ 
ferent spectral regions, i.e., in the CoRoT chromatic mode. 
The contamination factors (in a scale ranging from 0 to 1) 
are 0.026 and 0.036 in the two CoRoT runs, suggesting no 
relevant contribution from back- or fore-ground stars. 

The analysis of the original N2 data clearly detected 
a predominant peak at /i=0.3574 d“^ followed by several 
harmonics. Figure 1 shows LRaOl (top panel of the upper 
part) and LRa06 (top panel of the lower part) CoRoT N2 
data on a 15-d time baseline. All data points are plotted, 
including those acquired during the crossing of the South 
Atlantic Anomaly (SAA): they appear as bright outliers on 
a vertical line. We avoided any cleaning procedure since the 
detection of instrumental trends and jumps and their subse¬ 
quent correction is a crucial step to proceed in the analysis 
of the CoRoT data. Therefore, we computed a preliminary 
solution based on the set of frequencies /i, 2/i, 3/i, 4/i, 
and 5/i. The plots of the residuals (middle panels of the 
two parts of Fig. 1) show sudden jumps around JD 2454468 
in LRaOl and JD 2456003 in LRa06, the SAA outliers, and 
continuous oscillations. Outliers were iteratively removed at 
this stage. We still remained with the necessity to correct 
the sudden jumps and long-term drifts and to verify the 
origin of the residual oscillations. To do this, we applied 
the least-squares fits to each pulsation cycle, still using the 
set of frequencies from /i to 5/i. The 32-s cadence of the 
CoRoT data warranted more than 7000 measurements for 
each cycle, making the solution very robust. This procedure 
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Table 1. The Cepheids observed with CoRoT. CHR: chromatic mode, MON: white light only. 


CoRoT ID 

Other designation 

Run 

Start date 

End date 

AT 

Mode 

Cadence 

N 

0102618121 


LRaOl 

Oct. 23, 2007 

Mar. 3, 2008 

131.46 d 

CHR 

32 s 

326634 

0102618121 


LRa06 

Jan. 12, 2012 

Mar. 29, 2012 

76.61 d 

CHR 

32 s 

196032 

0300003751 

V376 Mon 

LRa02 

Nov. 16, 2008 

Mar. 8, 2009 

111.64 d 

CHR 

512 s 

16400 

0221644967 

DP Mon 

SRa02 

Oct. 11, 2008 

Nov. 12, 2008 

31.75 d 

CHR 

32 s 

72144 

0659466327 


LRc09 

Apr. 12, 2012 

July 5, 2012 

83.51 d 

MON 

512 s 

11894 

0315221415 


SRa03 

Mar. 5, 2010 

Mar. 29, 2010 

24.73 d 

MON 

512 s 

3670 

0211626074 


SRcOl 

Apr. 13, 2007 

May 9, 2007 

25.54 d 

MON 

512 s 

3996 

0221640090 

NSVS 12485452 

SRa02 

Oct. 11, 2008 

Nov. 12, 2008 

31.75 d 

CHR 

32 s 

73326 


returned the mean fluxes of the 47 cycles observed in the 
LRaOl and the 27 ones in LRa06 (Fig. 1, bottom panels of 
the two parts). These mean values show the variations due 
to the jumps and the long-term trend. We notice that the 
long-term trends behave in different ways in the two runs. 
The effect of the jump around JD 2454468 was noticeable 
in the data until JD 2454485, thus affecting the mean mag¬ 
nitudes of the subsequent six cycles. To remove the instru¬ 
mental variations, we interpolated the values of the mean 
fluxes (bottom panels of Fig.l) at the times of each individ¬ 
ual measurement and then we subtracted it. At this point 
two cleaned timeseries (326,634 and 196,032 points, respec¬ 
tively) were ready to be analysed in frequency. By applying 
this procedure we probably modified any possible additional 
signal around 2.8 d and this fact has to be carefully con¬ 
sidered when discussing the frequency content of the light 
curves of the two runs. 

We transformed the fluxes into magnitudes by calculat¬ 
ing the average flux of each run. We also measured the level 
of the noise of the CoRoT timeseries after subtracting the set 
of the /i frequency and harmonics. In the case of the LRaOl 
data we obtained 0.065 mmag in the 0.01-1.0 d“^ range, 
rapidly decreasing to 0.018 mmag in the 1.0-2.0 d~^ region, 
and more slowly to 0.012 mmag in the 2.1-3.0 d“^ region 
and to 0.007 mmag around 5.0 d“^. In the case of the LRa06 
data the four values are 0.104, 0.030, 0.020, and 0.010 mmag, 
respectively. The number of measurements (326,634 versus 
196,032) accounts for most of the difference between the 
noise levels of the LRaOl and LRa06 spectra. 
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Figure 2. Analysis of CoRoT N2 data of CoRoT 0102618121. 
Upper part, power spectra of the residuals after applying a global 
fit to the entire set of the LRaOl data (top panel) and LRa06 
data (bottom panel). Lower part: power spectra of the residuals 
after applying a fit to each pulsational cycle of the LRaOl data 
(top panel) and LRa06 data (bottom panel). 


4 A CRITICAL ANALYSIS OF THE 
FREQUENCY CONTENT OF 
CoRoT 0102618121 

The frequency analysis was performed by means of the iter¬ 
ative sine-wave, least-squares fitting method (Vanicek 1971) 
and then refined with the MTRAP algorithm (Carpino, Mi- 
lani & Nobili 1987), allowing us to include all the significant 
harmonics in the search for the best fit around /i. The best 
way to illustrate our final results is to follow the chronolog¬ 
ical order of their progressive consolidation. 

The analysis of the LRaOl data in white light was 
concluded in 2011 and the results were apparently very 
promising. We found the harmonics up to 6/i (ampli¬ 
tude 0.13 mmag), flanked by some modulation terms at 
±0.007 d~^ (i.e., 142 d). Since the separation is very close 
to the length of the observations, the origin of these terms 


has to be ascribed to a residual, still imperfect correc¬ 
tion of the instrumental drifts. Moreover, a few peaks not 
very close to /i and 2/i were detected: 0.488 d“^ (ampli¬ 
tude 0.40 mmag), 0.392 d“^ (0.14 mmag), and 0.288 d“^ 
(0.15 mmag). Only the peak at 0.488 d“^ had a S/N 
above the threshold considered as the significativity one in 
the Fourier analysis (S/N=4.0; Kuschnig et al. 1997), i.e., 
S/N=0.396/0.065=6.1. The power spectra illustrating the 
bunch of low amplitude peaks still centred around /i, the 
0.488-d“^ peak, and the 6/i harmonic are shown in the top 
panel of Fig. 2. However, the identification of the 0.488 d“^ 
peak as an additional mode was not really convincing be¬ 
cause no combination term was observed. 

Why this non-detection is so important? The combina¬ 
tion or coupling terms are generally ascribed to the non¬ 
linear responses of the stellar atmospheres to the strong ra¬ 
dial pulsations. The frequency analysis of the light curves of 
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Figure 3. Upper part, LRaOl data: light curves of the residuals 
after subtracting fi, 2/i, 3/i, 4/i, and 5/i calculated from the 
entire set of data (grey points) and from each pulsational cycle 
(black points). A part of the data is shown. Lower part, the same 
as in upper part, but for LRa06 data. 


double-mode Cepheids show several of them, also involving 
harmonics of the F and lO radial modes (Pardo & Poretti 
1997). The Fourier decomposition of the light curves also 
shows regularities in the distribution of the phase param¬ 
eters that are not yet fully understood (Poretti & Pardo 
1997). The combination terms have been detected in the 
light curves of other multimode pulsators, e.g. HADS stars 
(Poretti et al. 2005), RR Lyr stars (Gruberbauer et al. 2007), 
and white dwarfs (Wu 2001). The CoRoT photometry also 
allowed us to detect low-amplitude, nonradial modes in an 
HADS star pulsating in a single radial mode and also in this 
case a lot of combination terms were detected (Poretti et al. 
2011). Therefore, we should expect combination terms when 
additional modes, both radial and nonradial, are excited in 
a high-amplitude radial pulsator. 

The uncertainties on the real significativity of the peaks 
detected in the LRaOl data slowed down our theoretical 
analysis and made us unsure of how to proceed further. The 
re-observation of the same field in the LRa06 pointing sup¬ 
plied us with new data to verify the frequency content of 
CoRoT 0102618121. We applied the same procedure used 
for the LRaOl data to remove the content of the main vari¬ 
ation. The frequency analysis of the residuals of the LRa06 
data confirmed the presence of the 0.488 d“^ peak (Fig. 2, 
bottom panel of the upper part). Therefore, the confidence 
on the stellar origin of the 0.488 d“^ grew. However, the ra¬ 
tio 0.357/0.488=0.732 would be quite unusual for a double¬ 
mode Cepheid with a normal metallic content, both from ob¬ 
servational data (Pardo & Poretti 1997; Sziladi et al. 2007) 
and theoretical models (Buchler & Szabo 2007; Buchler 
2008). Therefore, such a ratio pointed out a possible nonra¬ 
dial mode. The identification of the peak as the rotational 


frequency is much more unlikely, since the corresponding ro¬ 
tational period (2.06 d) appears too short for a supergiant 
star. 

The excitation of nonradial modes in a classical Cepheid 
would be a very important ingredient in the modelling of 
the atmospheres of these supergiant stars, thus we submit¬ 
ted this result to further checks. At this stage, we detected 
the 0.488 d~^ peak in the residuals of each Long Run. These 
sets of residuals were obtained by applying a global fit to the 
LRaOl data (grey points in the upper panel of Fig. 3) and 
to the LRa06 ones (grey points in the lower panel of Fig. 3). 
Indeed, variations in the shape of the residuals are clearly 
visible, especially after the jumps at JD 2454468 and JD 
2456003. They remind the oscillations noticed before remov¬ 
ing the long-term trends and the jumps (Fig. 1, middle pan¬ 
els). This suggested that the computation of the mean fluxes 
was very effective in cleaning the power spectra, but the in¬ 
terplay betweeen the large amplitude of CoRoT 0102618121 
and the jump correction could somehow leave spurious vari¬ 
ability close to the significance level. For this reason we gave 
a closer look to the light curves of the residuals from which 
the 0.488 d“^ peak was detected. These light curves show 
the same type of oscillations, suggesting that the global so¬ 
lution was not able to match the instrumental alterations 
of the stellar signal. Thus we decided to calculate residuals 
following the same approach used to determine the mean 
fluxes. This implied to compute as many sets of parame¬ 
ters as the elapsed cycles were, i.e., 47 for LRaOl and 27 for 
LRa06. 

The analysis of these new sets of residuals were enlight¬ 
ening. The oscillations visible in Fig. 1 (middle panels) had 
completely disappeared and the light curves became really 
flat, with just some clues of high-frequency noise (Fig. 3, 
black points). Indeed, the power spectra of the residuals 
from fitting one cycle at a time show the peak at the 6/i 
harmonic only. This was observed both in the LRaOl data 
and in the LRa06 data (lower part of Fig. 2). Both spectra 
are characterized by an almost white noise. It could be ar¬ 
gued that the signal at 0.488 d“^ was removed when fitting 
the single cycles (also taking into account the almost exact 
3/4 ratio between the frequencies). However, the fact that 
the 6/i peak has remained clearly visible gave evidence that 
the fit with /i, 2/i, 3/i, 4/i, and 5/i of each cycle did not 
modified any signal component apart from these frequen¬ 
cies. We tried to come back to a new correction of the mean 
fluxes, but again the analysis of the global residuals thus 
obtained was unsatisfactory. We note that the technique of 
analysing single cycles were successfully applied to remove 
similar instrumental effects in the light curves of the high- 
amplitude 5 Set star CoRoT 101155310 (Poretti et al. 2011) 
and of the triple-mode Cepheid CoRoT 0223989566 (Poretti 
et al. 2014). 


5 THE PROPERTIES OF THE 

CoRoT 0102618121 LIGHT CURVES 

Once that we could rule out the excitation of other modes, 
we fitted the CoRoT magnitudes by means of a cosine se¬ 
ries (see appendix A). Table Al lists the coefficients of the 
solutions of the light curves of both runs. The shape of the 
light curve could be defined by the Fourier parameters (see 
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Table 2. Fourier parameters of the curves of CoRoT 0102618121 in the different colours. 


Term LRaOl LRa06 



white 

rc 

ac 

be 

white 

rc 

ac 

be 

021 [rad] 

4.1101(7) 

4.220(2) 

4.070(2) 

3.934(2) 

4.1081(9) 

4.191(2) 

4.012(6) 

3.925(6) 

031 [rad] 

4.645(6) 

4.637(13) 

4.499(21) 

4.683(20) 

4.689(7) 

4.674(15) 

4.814(84) 

4.920(64) 

R 21 

0.05395(4) 

0.05371(9) 

0.05260(12) 

0.05166(11) 

0.05646(5) 

0.05717(14) 

0.05265(34) 

0.04961(30) 

041 [rad] 

3.624(7) 

3.857(18) 

3.431(21) 

3.240(17) 

3.573(7) 

3.712(26) 

3.290(34) 

3.193(32) 

Full ampl. [mag] 

0.222 

0.191 

0.250 

0.306 

0.233 

0.208 

0.288 

0.364 

Res. rms [mag] 

0.0017 

0.0036 

0.0061 

0.0063 

0.0018 

0.0045 

0.0150 

0.0165 


CoRoT 0102618121 



4400 4450 4500 5950 6000 

HJD-2450000 


Figure 4. The Fourier parameters and O-C values as measure¬ 
ments of the cycle-to-cycle variations in the CoRoT 0102618121 
light curve. The grey boxes are the ibScr bands around the average 
values. Left part: LRaOl data; Right part: LRa06 data. 

appendix A) and Table 2 lists those relevant to our discus¬ 
sion. 

5.1 The stability of the light curves 

The absence of any relevant peak once that the contribution 
of the main oscillation has been removed is not by itself an 
evidence that the light curve is stable. Cycle-to-cycle varia¬ 
tions as those observed in the Kepler light curve of VI154 
Cyg or in the radial velocity curves could not result in a well 
defined peak in the power spectrum. Moreover, at least in 
principle, the residual peaks left around fi when applying 
a global fit could be ascribed to cycle-to-cycle variations, 
maybe erratic ones. 

We analysed the R 21 , 02 i, and (j> 3 i parameters of each 
cycle to quantify how much the light curve deviates from the 
average behaviour. The CoRoT 32-s cadence allowed us to 
get about 7500 measurements per cycle, adequate to survey 
the harmonics 2/i and 3/i in a very accurate way, though 
their amplitudes are very small (6.0 and 0.7 mmag, respec¬ 
tively). The Fourier parameters are shown in Fig. 4. The 
error bars are three times the standard deviations (±3 (j) 
around the average of the parameters of the individual cy¬ 
cles. The immediate result is that the differences from the 
average values are small. We could notice that all the indi¬ 


vidual values are between 3cr, though some apparently regu¬ 
lar variations are sometimes visible (e.g. the R 21 behaviour 
in the LRaOl). Therefore, we performed a frequency analysis 
of the Fourier parameters. We could not detect any period¬ 
icity since the power spectra of the R 21 , 021 , and 03 i pa¬ 
rameters do not show any common feature either in a single 
Long Run or when comparing LRaOl and LRa06 results. 

However, the best parameter to quantify the sta¬ 
bility of the light curve is probably the O-C parame¬ 
ter. We determined the observed Tmax’s by means of cu¬ 
bic spline functions with a non-zero smoothing param¬ 
eter that depends on the number of measurements to 
fit and on their scatter (Reinsch 1967). We thus ob¬ 
tained the linear ephemeris from the least-squares method 
Tmar, = HJD 2454400.2977 -b 2.798061 • E 

±0.0008 ± 0.000002 

(E= 0,1, 2,..., n) and hence we could calculate the O-C value 
for each Tmax- The residual rms of the least-squares linear 
fit is 0.0052 d and the 3a band is shown in Fig. 4. We also 
estimated the uncertainty on the individual Tmax determi¬ 
nation as the difference between the two times correspond¬ 
ing to the intersection of the spline function with the line 
y = rrimax + ojyjn — 1 (Fig. 5), where mm.ax is the mag¬ 
nitude at maximum and cr is the standard deviation of the 
fit. Such uncertainties are ranging from 0.010 to 0.016 d and 
hence the 3a interval (i.e., 22.5 min) matches the individual 
error bars, too. We analysed in frequency the O-C values 
too, again obtaining no clear indication of periodicities. 

As a final conclusion, looking at the different plots 
shown in Fig. 4 and at the results of the various frequency 
analyses, we cannot identify any additional periodicity in the 
variability of CoRoT 0102618121. The Fourier parameters 
and the O-C values sometimes mimic an apparent regular 
behaviour. This could be a hint of cycle-to-cycle variations, 
much less regular than those observed in the Kepler light 
curve of V1154 Cyg (Kanev et al. 2015). 


5.2 The comparison between LRaOl and LRa06 
light curves 

We merged the LRaOl and LRa06 measurements and we 
calculated a refined value of the period on the whole 
time baseline. The MTRAP code supplied the final vale 
P=2.7980596±0.0000001 d, in excellent agreement with the 
value obtained from the Tmax’s values only. The most rel¬ 
evant difference in the datasets of the two Long Runs is a 
slight increase from 0.222 to 0.233 mag in the full-amplitude 
of the light curves. Since we are dealing with absolute pho- 
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HJD-2454399 

Figure 5. Determination of the error bar ATmax of the first 
Tmax of CoRoT 0102618121. The original data are in grey, the cu¬ 
bic spline fit is the blue continuous line, the y = nimax+crlViT- ~ 1 
the horizontal dashed line; the two red circles indicate the inter¬ 
section between the spline curve and the line (better shown in the 
inset). The differences of their abscissae with respect to the Tmax 
one fix the error bars ATmax. 

tometry, i.e., without any comparison star, we have to be 
sure that the measurements were performed in the same in¬ 
strumental system to infer that the amplitude difference has 
a stellar origin. We immediately noticed that the star was 
observed in the same CCD in the two runs, but with a differ¬ 
ent mask and at different pixel coordinates. As a results, the 
average white flux measured was of 1 030 069 electrons in 
LRaOl and 944,911 electrons in LRa06. Moreover, the flux 
subdvision in rc, gc, be colours was also different: 65, 14, 21 
per cent in LRaOl, 78, 10, 12 per cent in LRa06. We conclude 
that the small difference in amplitude could very probably 
be ascribed to the different settings of the CoRoT pointings 
in the two runs. The same effect occurs in Kepler photome¬ 
try: average fluxes and amplitudes change for the same star 
observed in the same CCD but with different masks. These 
effects are particularly disturbing when dealing with large 
amplitude variables (Benko et al. 2014). 

The procedure applied to white light measurements was 
repeated for each colour of both runs. Other similar, sub¬ 
tle instrumental effects became evident during this process 
(Poretti, Le Borgne, & Klotz 2015). Then, the Fourier de¬ 
composition was applied to the cleaned timeseries in each 
colour (Table 2). The parameters clearly show the expected 
dependence from the bandpass (see Morgan, Simet & Bar- 
genquast 1998, for a comparison between V and I param¬ 
eters). The different flux distributions between LRaOl and 
LRa06 also account for the small differences between param¬ 
eters in the same colour. 


6 DATA ANALYSIS OF THE OTHER 
CEPHEIDS OBSERVED WITH COROT 

The procedure described in the previous sections for 
CoRoT 0102618121 was applied to the other Cepheids, 
whose folded light curves are shown in Fig. 6 (six Pop. I 


stars) and Fig. 7 (one Pop. II star). The Fourier parameters 
are listed in Table 3. 


6.1 Stars observed in the Long Runs 

The frequency analysis of the V376 Mon data supplied 
/i=0.6060 d"\ refined to P = 1.651969 ± 0.000001 d with 
MTRAP. Table Al lists the least-squares solution of the 
CoRoT white light curve. The frequency analysis did not 
detect any significant additional mode, but again the resid¬ 
ual power spectra strongly depend on the procedure used to 
prewhiten the data. We obtained a bunch of peaks around 
fi when applying a global fit up to 7/i to the entire dataset. 
These peaks disappeared when the fit was calculated cycle- 
by-cycle (Fig. 8, upper part). The harmonic 8/i was de¬ 
tected in both cases, confirming that the single cycle fit 
did not alter the frequency content. The Fourier parame¬ 
ters and the O-C values show a small scatter around the 
means (Fig. 9), suggesting no periodic variation. 

The frequency analysis of CoRoT 0659466327 data re¬ 
turned /i=0.648 d“^ and the period was refined to P = 
1.542353 ± 0.000002 d. Table A2 lists the least-squares so¬ 
lution of the white light curve. The bump on the ascending 
branch (Fig. 7) is typical of Pop. II Cepheids. Differently 
than in the previous cases, the peak in the power spec¬ 
trum of the residuals from the global fit up to 9/i occurred 
around 2/i, but the bunch of peaks at /i were detected, too. 
The power spectrum of the residuals from the cycle-by-cycle 
fits is completely flat and it just shows the lO/i harmonic, 
not prewhitened from the data (Fig. 8, lower part). Small- 
size scale variations are observed in the Fourier parameters 
and O-C values, but again without any trace of periodicity 
(Fig. 10). 

The power spectra of the residuals of 
CoRoT 0102618121, V376 Mon and CoRoT 0659466327 
calculated by means of a global least-squares solution show 
the limits of this approach when applied to continuous and 
high-precision timeseries where systematic errors (jumps 
and/or drifts) can occur. When possible, e.g. in presence of 
a high-amplitude predominant mode, the frequency analysis 
of cycle-to-cycle residuals are preferable, taking care to 
keep a known frequency, e.g. an harmonic, as internal check 
of the signal preservation. 


6.2 Stars observed in the Short Runs 

The solutions of the light curves of the Cepheids 
CoRoT 0221640090 (Table A3), CoRoT 0211626074 (Ta¬ 
ble A4, CoRoT 0315221415 (Table A5), and DP Mon (Ta¬ 
ble A6 were much less laborious since these stars were ob¬ 
served once in some Short Runs (SRs) and hence followed 
for a restricted number of cycles. The O-C plots do not pro¬ 
vide useful hints about periodic variations (Fig. 11). Only 
in the case of CoRoT 0221640090 (bottom panel) the string 
of consecutive Tmax shows a hint of regular behaviour, but 
the monitoring is too short to assess its repetitivity on a 
long time-scale. Moreover, the standard deviation around 
the linear ephemeris is ±0.0033 d and all the maxima are 
well within 3cr. 
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Figure 6. Family portrait of the CoRoT light curves of Pop. I Cepheids. 


Table 3. Fourier parameters of the curves of the Cepheids observed with CoRoT. See Table 2 for CoRoT 0102618121. 


Star 

Mode 

Tmax[HJD- 

2450000] 

Period 

[d] 

Ampl. 

[mag] 

021 

[rad] 

R 21 

031 

[rad] 

</>41 

[rad] 

V376 Mon 

F 

4785.4731 

1.651969(1) 

0.372 

4.1811(3) 

0.3745(1) 

2.2206 (1) 

0.172(1) 

CoRoT 0221640090 

10 

4754.6755 

2.19758(4) 

0.288 

4.560(1) 

0.1397(1) 

3.641(5) 

2.670(9) 

DP Mon 

F 

4754.0407 

3.08704(4) 

0.678 

4.312(1) 

0.5201(1) 

2.4703(4) 

0.5439(6) 

CoRoT 0315221415 

F 

5258.4566 

4.7085(17) 

0.557 

4.363(1) 

0.4650(5) 

2.478(3) 

0.434(5) 

CoRoT 0211626074 

F 

4199.7103 

5.4706(2) 

0.296 

4.731(1) 

0.3147(3) 

3.242(3) 

1.391(35) 

CoRoT 0659466327 

Pop. II 

6031.4482 

1.542353(2) 

0.492 

4.394(1) 

0.2913(2) 

0.886(1) 

4.953(2) 


7 THE FOURIER PARAMETERS AS 
PULSATION MODE DISCRIMINANT 

The Fourier decomposition of the light curves of Cepheids al¬ 
lowed us to disentangle F and lO radial pulsators (Fig. 12). 
The (() 2 i-sequence formed by the classical Cepheids repro¬ 
duces the well-known Hertzsprung progression. It can be 
described by the linear relation <(>21 = 3.332 -|- 0.216 P 
(Poretti 2000). The progression is sharp, since a shift of 
±0.30 rad pushes a star well outside it. Another Z-shaped 
sequence was put in evidence by applying the Fourier de¬ 
composition to the light curves of Galactic Cepheids with 
P <8 d: the upper part of this sequence is composed of stars 


with 021 >4.2 rad, the lower part of stars with 02i< 4.0 rad 
(Fig. 12, top panel). The (f> 3 i—P plane was the Rosetta stone 
to link the upper and lower sequences, since it shows just one 
continuous sequence, well separated from the classical one 
(Fig. 12, middle panel). The Fourier decomposition of the 
light curves of double-mode (F/IO) Cepheids provided the 
definitive evidence that the pulsation mode drives the sep¬ 
aration between the different sequences (Pardo & Poretti 
1997; Poretti & Pardo 1997). The R 21 —P plot shows that 
the amplitude of the harmonic 2/i of the lO-mode pul¬ 
sators is much smaller than that of the main frequency fi 
(Fig. 12, bottom panel). This fact makes the light curves of 
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Figure 7. The light curve of the Pop. II Cepheid CoRoT 
0659466327. 
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Figure 8. Upper part-, analysis of CoRoT N2 data of V376 Mon: 
power spectra of the residuals after subtracting a global fit up 
to 7/i from the entire set of data (top panel), and a single fit 
from each pulsational cycle (bottom panel). Lower part-, the same 
analysis for CoRoT 0659466327, with a fit up to 9/i. 


the lO-mode pulsators much more (but not exactly) sine¬ 
shaped than those of F-mode pulsators. Theoretical works 
demonstrated how the inclusion of convective energy trans¬ 
port has made possible to reproduce the observed sequence 
of lO-pulsators and that the sharp features in the Fourier 
coefficients are due to the resonance Pi /P4 = 2 between the 
fourth and first overtones (Feuchtinger, Buchler & Kollath 
2000 ). 

The classification of CoRoT 0211626074 and 
CoRoT 0315221415 is straightforward, since the (j) 2 i, 
R 21 , and (f> 3 i parameters put both these long-period (in our 


CoRoT 0300003751-V376 Mon 



Figure 9. The Fourier parameters and O-C values as measure¬ 
ments of the cycle-to-cycle variations in the V376 Mon light 
curve. The grey boxes are the ±30- bands around the average 
values. 
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HJD-2450000 


Figure 10. The Fourier parameters and O-C values as measure¬ 
ments of the cycle-to-cycle variations in the CoRoT 0659466327 
light curve. The grey boxes are the ±3fT bands around the average 
values. 


sample) Cepheids on the well-defined sequence of P-mode 
pulsators. 

The period of CoRoT 0102618121 and DP Mon are very 
particular in this context. They lie exactly on the disconti¬ 
nuity of the Z-shaped sequence of the lO pulsators in the 
4>2i—P plane, where the intersection with the P-sequence 
occurs. Therefore, the (l> 2 i parameter alone is not able to 
identify their pulsation mode. However, the much higher 
value of <j) 3 i and the much smaller value of P 21 with respect 
to those of the classical Cepheids clearly demonstrate that 
CoRoT 0102618121 is a lO pulsator. On the other hand, the 
(j) 3 i and R 21 parameters of DP Mon match the P-sequence. 
The fight curve of DP Mon is very asymmetrical (Fig. 6) 
and it is very similar to those of RRab stars. In particular, 


MNRAS 000, 1-12 (2015) 

















































10 E. Poretti et al. 


"d 

u 

o 



1 ^ 

T 



1 

-1 1 1 1 t- 

CoRoT 0315221415 

7 

- 

1 • 




• 


- 

: 

_1_L 

_1 



J_ 

j_^^_1_ 

J 


5260 5270 52B0 5290 5300 


u 


o 


^ ^ r 

^ ^ ^ 



CoRoT 0211626074 C 



- 


• • 

- 

-_^_L 

_^_1_,_ 1 _ ^ _ 

_^_1_^_- 


4200 4210 4220 4230 4240 



4750 4760 4770 47B0 4790 

HJD-2450000 


Figure 11. The O-C values with respect to a linear 
ephemeris fitting the Tmax of CoRoT 0221640090, DP Mon, 
CoRoT 0315221415, and CoRoT 0211626074. The grey box is the 
ibSfT band (where cr is the standard deviation of the fit) around 
the average values. 



Log Period 


Figure 12. The separation between fundamental (open green 
circles) and first overtone (filled black circles) Galactic Pop. I 
Cepheids in the progression of the Fourier parameters (in V light) 
in function of the period. The red squares indicate the posi¬ 
tions of the CoRoT Cepheids: V376 Mon, CoRoT 0221640090, 
CoRoT 0102618121, DP Mon, CoRoT 0315221415, and 
CoRoT 0211626074 (from left to right). 


Table 4. Phase shifts of the light curves of the Cepheids observed 
in the CoRoT chromatic mode. 


Star 

Blue 

Green 

Red 

White 

CoRoT 0102618121 

0.993 

0.997 

1.004 

1.000 

in LRaOl 

±0.003 

±0.003 

±0.003 

±0.003 

CoRoT 0102618121 

0.994 

0.996 

1.003 

1.000 

in LRa06 

±0.005 

±0.005 

±0.003 

±0.003 

CoRoT 0221640090 

0.994 

1.002 

1.007 

1.000 


±0.003 

±0.003 

±0.005 

±0.003 

V376 Mon 

0.996 

0.998 

1.002 

1.000 


±0.003 

±0.003 

±0.003 

±0.003 

DP Mon 

0.998 

1.000 

1.001 

1.000 


±0.003 

±0.003 

±0.003 

±0.003 


notice the kink before the minimum, that requires a large 
number of harmonics to be fitted. 

We classify CoRoT 0221640090 as a lO pulsator, like 
CoRoT 0102618121. The case of V376 Mon is more intrigu¬ 
ing, since its period is very short also with respect to the full 
sample of Galactic Cepheids. The Cepheids with P <2.5 d 
are quite rare and the space of Fourier parameters is scarcely 
populated in this region. For this reason the values of 4>2i 
and (/) 3 i are not able to provide a clear mode identification. 
However, the R 21 value is quite high, and hence the star is 
very probably a T—mode pulsator. 


T.l The coloured light curves 

The continuous photometry allowed us to verify the phases 
of maximum brightness (0max) in the CoRoT colours with 
great accuracy. The analyses described in the previous sec¬ 
tions were repeated for each colour of each star. The results 
are summarized in Table 4. There is a common trend indi¬ 
cating that the maximum brightness occurs a little earlier 
at short wavelengths than at long ones. The only exception 
is DP Mon, whose maxima are measured at the same phases 
(within error bars) in all colours. 

The observed phase shifts are very small, around 
0.01 period. This result strongly corroborates the practice 
to merge all the maxima observed with different filters or 
techniques when studying secular period variations. Indeed, 
this shift is usually smaller than the uncertainties on the 
Tmax determination in case of ground-based, fragmented ob¬ 
servations. This is the first measurement of the phase shifts 
among light curves of Cepheids in different colours after the 
pioneering approach by Wisniewski & Johnson (1968). 


8 DISCUSSION AND CONCLUSIONS 

The power spectra of the CoRoT light curves of Galactic 
Gepheids did not show any convincing evidence of the excita¬ 
tion of nonradial modes, especially in the three monitored in 
the Long Runs. Stars showing additional modes were found 
with an incidence of 10% among the lO-Cepheids of the 
LMC (Moskalik & Kolaczkowski 2008). The only two lO- 
pulsators observed by GoRoT cannot yet tell us if Galac¬ 
tic lO-Gepheids behave differently than LMG ones. On the 
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other hand, the F-Cepheids of the LMC do not show ad¬ 
ditional modes (Moskalik & Kolaczkowski 2008, 2009) and 
the intensive CoRoT photometry does not modify this is¬ 
sue. The detection of the same peak in two independent sets 
of residuals of CoRoT 0102618121 seemed to be a decisive 
piece of evidence in favour of the excitation of a nonradial 
mode, but it not survived the change of the procedure for 
the computation of the residuals. 

In our opinion the instrumental effects seem to play an 
important role in the long time series of space photometry. 
Our recommendation is to analyse residuals calculated from 
local solutions (e.g. cycle-to-cycle when possible) to search 
for additional periodicities. Instrumental drifts and jumps 
(Fig. 1) or, in extreme cases, numerical errors or amplitude 
changes can enhance spurious peaks in the power spectrum 
of the residuals calculated from a global fit (Fig. 2). The 
residuals from the cycle-by-cycle fits are providing a much 
more reliable dataset to be analysed in frequency (Fig. 3). 
Omitting a high harmonic from the fit and verify that it is 
correctly recovered is an excellent litmus test of the whole 
procedure. Moreover, the lack of combination terms between 
the main pulsation mode and the suspected additional com¬ 
ponent should be considered a warning against the identifi¬ 
cation of this component as a true excited mode. 

Another compelling result expected from the analysis of 
the CoRoT light curves was to verify the repetitivity and the 
amplitude of the cycle-to-cycle variations, as those observed 
in the Kepler data of V1154 Cyg (Kanev et al. 2015) and in 
the radial-velocity surveys (Anderson 2014). The Fourier pa¬ 
rameters of the light curves of CoRoT Cepheids show some 
variability, but in general all the changes are within ±3 ct. 
Moreover, they seem to be erratic, not reproducing phenom¬ 
ena like activity cycle or long-term effects. 

No large scatter was observed around the maxima 
of the folded light curves and no large variations were 
observed in the O-C plots of the Cepheids monitored 
with CoRoT. In particular, the O-C values are all within 
±0.01 d, i.e.,±15 min, including the lO-mode pulsators 
CoRoT 0102618121 and CoRoT 0221640090 (Figs. 4, 9, 10, 
and 11). This is particularly interesting since the space pho¬ 
tometry obtained with the satellite MOST showed some in¬ 
stabilities in the pulsation cycle of the lO-Cepheid SZ Tan, 
while the F-mode RT Aur repeated the light curve more pre¬ 
cisely (Evans et al. 2015). Therefore, we performed some ad¬ 
ditional tests on the folded light curves, evaluating the stan¬ 
dard deviations as a function of phase bins. Tiny increases 
of the standard deviations were measured in correspondence 
of the most rapid light variation (i.e., on the steep ascend¬ 
ing branch); they mostly reflect the variability due to the 
pulsation within the phase bins, as proven by the fact that 
they appear to be more relevant for the stars showing the 
largest amplitudes and the most asymmetrical light curves. 
This suggests that the standard deviations of the phase bins 
on the rising branch cannot be used as a tool to search for 
cycle-to-cycle variations. We can conclude that in the case 
of the seven CoRoT Cepheids, the cycle-to-cycle variability 
is very limited and not greatly affecting the repetitivity of 
the light curves. The light variations observed around the 
maxima of SZ Tan are more promising, but on this context 
the low-frequency peaks in the power spectrum have to be 
carefully considered: the residuals observed in MOST pho¬ 


tometry (Fig. 6 in Evans et al. 2015) are very similar to the 
spurious ones detected in CoRoT one (Fig. 3). 

We also benefit from the unique opportunity provided 
by the CoRoT multicolour photometry. Very small shifts in 
the Tniax’s were observed for the first time in the optical 
range. Because these shifts are around 0.01 period only, our 
analysis support the use of Tmax’s observed in different pass- 
bands to reconstruct secular period variations (e.g. Evans 
et al. 2015). 

The accurate values of the Fourier parameters provided 
by the available space photometry up to very high harmon¬ 
ics can be used to improve the automatic classification of the 
Cepheid subclasses, i.e., Pop. I, Pop. II, F-mode, and lO- 
mode pulsators. The quantitative separation of these sub¬ 
classes in the parameter space is of particular relevance when 
considering the large harvest of new variables that will be 
discovered by future, wide-held, and all-sky space missions 
like TESS (Ricker et al. 2014) and PLATO 2.0 (Rauer et al. 
2014). 
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APPENDIX A: SOLUTIONS OF THE COROT 
LIGHT CURVES 

We provide the solutions of the CoRoT light curves of the 
Cepheids observed in the Long Runs. We fitted the CoRoT 
magnitudes by means of the formula 

m{t) = mo + Y cos[27r i f{t — To) + cf)i] (Al) 


and we calculated the Fourier parameters Rij = AijAj (e.g. 
R 21 = A 2 /A 1 ) and (j>ij = j (j>i- i (j)j (e.g. <^21 = <j >2 - 2 (pi 
and <^31 = (j )3 — 3 (pi) from the amplitude {Ai) and phase 
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Table Al. Coefficients and formal errors (Ifj) of the least—squares solutions of the light curves of the Pop. I, lO-mode pulsator 
CoRoT 0102618121 and of the Pop. I, i^-mode pulsator V376 Mon. 


Term 

CoRoT 0102618121, /i =0.3573905(1) d"! 

V376 Mon 

/i=0.6053382(4) d"! 

LRaOl 

LRa06 

Ampl. 

[mag] 

Phase 

[rad] 

Ampl. 

[mag] 

Phase 

[rad] 

Ampl. 

[mag] 

Phase 

[rad] 

h 

0.110500(4) 

3.06375(4) 

0.116054(6) 

3.05962(5) 

0.15852(2) 

2.5446(1) 

2/l 

0.005962^) 

3.9544o52) 

0.0065525) 

3.944145o) 

0.059375) 

2.987i5) 

3/i 

0.0006995) 

1.269655i6) 

0.0008965) 

1.30168555) 

0.027155) 

3.57125) 

4/l 

0.00062i5) 

3.31223(693) 

0.0008095) 

3 . 24492526 ) 

0.010965) 

4.067654) 

5/i 

0.0002875) 

4.55196(1499) 

0.0003895) 

4 . 519495509 ) 

0.004355) 

4 . 450354 ) 

6/i 

0.0001355) 

5.74428(3180) 

0.00012o5) 

5.64753 5876) 

0.001565) 

4 . 668954 ) 

7/i 

0.0000495) 

0.55276 5822) 

0.0000785) 

0.098525486) 

0.000775) 

4 . 5913590 ) 

8/i 

0.000024(4) 

2 . 1097957710 ) 

0.0000255) 

1.9764753114) 

0.000495) 

4.6907597) 

9/i 





0.00036(2) 

4.92985o8) 

lO/i 





0.000285) 

5.357o534) 

ll/i 





0.000165) 

5 . 86 I 5552 ) 

12/i 





o.oooio5) 

6.16335511) 

13/i 





0.000075) 

0.759o5i55) 

14/i 





0.000055) 

0.74225o75) 

Res. rms 

0.00174 mag 

0.00184 mag 

0.00133 mag 

To 

HJD 2454394.6990 

HJD 2455936.4282 

HJD 2454785.4731 


Table A2. Coefficients and formal errors (Icr) of the least- 
squares solution of the light curve of the Pop. II Cepheid 
CoRoT 0659466327. 



CoRoT 0659466327 
/i =0.648366006 d"! 

Term 

Ampl. 

Phase 


[mag] 

[rad] 

fl 

0.19348(4) 

2.6323(2) 

2/1 

0.056375) 

3.37525) 

3/1 

0.034785) 

2.5001(12) 

4/1 

0.027055) 

2.915556) 

5/1 

0.016085) 

3.8907(26) 

6/1 

0.011545) 

4.795i56) 

7/1 

0.006295) 

5.6I7357) 

8/1 

0.002625) 

0.3086(161) 

9/1 

0.000555) 

1.1627567) 

10/1 

0.000525) 

6.0158532) 

11/1 

0.000765) 

0.8167555) 

12/1 

0.000835) 

1.9510(508) 

13/1 

0.000685) 

2.95675i4) 

14/i 

0.0004i5) 

4.14195o24) 

15/1 

0.000195) 

5.20845i96) 

I 6/1 

0.000075) 

5.7712(6216) 

17/1 

0.000095) 

5.42645543) 

I 8/1 

0.000165) 

0.26525712) 


Res. rms 0.00323 mag 

To HJD 2456031.4482 


Table A3. Coefficients and formal errors (Icr) of the least- 
squares solution of the light curve of the Pop. I Cepheid, lO-mode 
pulsator CoRoT 0221640090. 



CoRoT 0221640090 


/i=0.455047 d-i 

Term 

Ampl. 

Phase 


[mag] 

[rad] 

fl 

0.14079(1) 

2.9111(1) 

2/1 

0.019675) 

4.09875) 

3/1 

0.0029o5) 

6.091o57) 

4/1 

0.001565) 

1.7476 57 ) 

5/1 

0.001025) 

3.1642532) 

6/1 

0.000555) 

4.3818544) 

7/1 

0.000355) 

5.517o584) 

8/1 

0.000175) 

0 . 4274595 ) 

9/1 

0.000095) 

1 . 39755468 ) 

10/1 

0.000025) 

1.96385i12) 

11/1 

0.000065) 

3.291o5i57) 

12/1 

0.000025) 

3.62095i19) 

Res. rms 

0.00259 mag 

To 

HJD 2454754.6755 
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Table A4. Coefficients and formal errors (Icr) of the least- Table A6. Coefficients and formal errors (Icr) of the least- 

squares solution of the light curve of the Pop. I Cepheid, F-mode squares solution of the light curve of the Pop. I Cepheid, F-mode 

pulsator CoRoT 0211626074. pulsator DP Mon. 



CoRoT 0211626074 
/i=0.182796 d-i 

Term 

DP Mon 
/i=0.323939 d-i 

Term 

Ampl. 

[mag] 

Phase 

[rad] 

Ampl. 

[mag] 

Phase 

[rad] 

fl 

0.12984(3) 

2.5207(3) 

fl 

0.23679(2) 

2.2927(1) 

2/l 

0.04086(3) 

3.48975) 

2/1 

0.123155) 

2.61455) 

3/l 

0.01028(3) 

4.5213(32) 

3/1 

0.071985) 

3.06525) 

4/l 

0.00095(3) 

5.1938 541) 

4/1 

0.040445) 

3.43155) 

5/l 

0.001665) 

3.0263595) 

5/1 

0.023525) 

3.69965) 

6/i 

0.001285) 

3.7305553) 

6/1 

0.014445) 

4.003652) 

7/i 

0.001035) 

4.58255i5) 

7/1 

0.008565) 

4.271i5o) 

8/i 

0.000625) 

5 . 3574527 ) 

8/1 

0.005235) 

4.4842 53) 

9/i 

0.000365) 

6.26635i2) 

9/1 

0.003235) 

4 . 687254 ) 

lO/i 

0.000185) 

0.8783 5842) 

10/1 

0.002055) 

4.789755) 

ll/i 

o.oooio5) 

2.12665298) 

11/1 

0.001265) 

4 . 8674537 ) 

12/i 

0.000085) 

2.79985i24) 

12/1 

0.000965) 

5.0046580) 




13/1 

0.000745) 

5 . 0513534 ) 

Res. rms 

U.UU144 mag 

14/1 

0.000585) 

5 . 2112596 ) 

To 

HJD 2454199.7103 

15/1 

0.000485) 

5.3839559) 




16/1 

0.0004i5) 

5.632i525) 




17/1 

0.000345) 

5.975i5i5) 




18/1 

0.000245) 

6.08235i2) 

Table A5. 

Coefficients 

and formal errors (Icr) of the least- 

19/1 

0.000225) 

0.1627599) 

squares solution of the light curve of the Pop. I Cepheid, F-mode 

20/1 

0.000145) 

0.52565234) 

pulsator CoRoT 0315221415. 

21/1 

0.000165) 

0.72885i13) 




22/1 

o.oooio5) 

1.360i5793) 


CoRoT 0315221415 

23/1 

0.000095) 

1.187i5o08) 


fi=0.212384 d-i 

24/1 

0.000035) 

1 . 47775286 ) 

Term 

Ampl. 

Phase 

25/1 

0.000055) 

1.81885629) 


[mag] 

[rad] 

Res. rms 

0.00328 mag 

fl 

0.21116(6) 

2.4327(3) 

To 

HJD 2454754.0407 

2/l 

0.098195) 

2.94565) 




3/i 

0.044155) 

3 . 492754 ) 




4/l 

0.0201i5) 

3.88195i) 




5/i 

0.007355) 

4.091654) 




6/i 

0.003865) 

3.9842559) 




7/i 

0.002985) 

3.70865o7) 




8/i 

0.003145) 

3.896i596) 




9/i 

0.002995) 

4.3491(205) 




lO/i 

0.002465) 

4.8484548) 




ll/i 

0.0020i5) 

5.42955o5) 




12/i 

0.001585) 

6.0745 5 87) 




13/i 

0.001225) 

0.37695o1) 




14/i 

0.000875) 

I.I 6 I 25 O 6 ) 




15/i 

0.0006o5) 

1.81275o16) 




16/i 

0.000455) 

2.672o5363) 




17/1 

0.000325) 

3.5748 5915) 




18/i 

0.000165) 

4.250o5960) 




19/i 

0.0002o5) 

5.16835o62) 




20/i 

0.000165) 

0.3755 5767) 





Res. rms 0.00262 mag 

To HJD 2455258.4566 


MNRAS 000, 1-12 (2015) 

















